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Abstract 

We discuss various issues associated with the calculation of the reduced functional determi- 
nant of a special second order differential operator F — ~S jdrr^ + gjg, g = d^g/dr^, with a 
generic function (?(t), subject to periodic and Dirichlet boundary conditions. These issues include 
the gauge- fixed path integral representation of this determinant, the monodromy method of its 
calculation and the combination of the heat kernel and zeta-function technique for the derivation 
of its period dependence. Motivations for this particular problem, coming from applications in 
quantum cosmology, are also briefly discussed. They include the problem of microcanonical ini- 
tial conditions in cosmology driven by a conformal field theory, cosmological constant and cosmic 
microwave background problems. 

1. Introduction 

Several essays on calculational methods in quantum problems, we are honored to present here, are ded- 
icated to Professor Dowker to praise his fundamental and thought-provoking contributions to quantum 
field theory in curved spacetime, spectral geometry and various methods of mathematical physics. His 
results in quantum theory with external gravitational and matter fields, pioneering contribution to the 
calculational method of zeta-function, treatment of zero modes, conformal anomalies and boundary 
terms, etc. determine the scope of methods and issues that we are going to discuss here. 

The importance of these methods follows from a simple fact that successful applications on non- 
trivial backgrounds always present a challenge and can be accomplished in a closed analytic form 
only in a limited number of cases. This is equally true with regard to calculations in the quantum 
mechanical sector of field models even despite simplifications occurring in this sector due to its spatial 
homogeneity or other symmetries. A particular case of such calculations is the class of problems 
involving the differential operator of the form 



(1-1) 

where g = gir) is a rather generic function of its variable r. From calculational viewpoint, the virtue 
of this operator is that g{T) represents its explicit basis function - the solution of the homogeneous 
equation, 

Fg{T) = 0, (1.2) 
which immediately allows one to construct its second linearly independent solution 

*M = »m/^ (1.3) 

r * 

and explicitly build the Green's function of F with appropriate boundary conditions. On the other 
hand, from physical viewpoint this operator is interesting because it describes long-wavelength per- 
turbations in early Universe, including the formation of observable CMB spectra [T] [2], statistical 
ensembles in quantum cosmology [3], etc. In particular, for super horizon cosmological perturbations 
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of small momenta fc^ ^ g/g their evolution operator only slightly differs from (|l.ip by adding fc^ 
to its potential term, whereas in the minisuperspace sector of cosmology, corresponding to spatially 
constant variables, the operator has exactly the above form. 

Up to an overall sign, this operator is the same both in the Lorentzian and Euclidean signature 
spacetimes with the time variables related by the Wick rotation r — it. In the Euclidean case it 
plays a very important role in the calculation of the statistical sum for the microcanonical ensemble 
in cosmology. This ensemble realizes the concept of cosmological initial conditions by generalizing the 
notion of the no-boundary wavefunction of the Universe |4] to the level of a special quasi-thermal state 
and puts it on the basis of a consistent canonical quantization [5j |6] . This concept is very promising 
both from the viewpoint of foundations of quantum cosmology and their applications within the 
cosmological constant, inflation and dark energy problems ^ IHl [3 [HI [3] ■ 

The recently suggested alteration in foundations of quantum cosmology - the theory of initial 
conditions for the early Universe - consists in a qualitative extension of the class of its initial quantum 
states. Instead of a usually accepted pure state, like the no-boundary one, it is assumed that the 
cosmological state can be mixed and characterized by the density matrix [5]. Under a natural and 
most democratic assumption of the microcanonical distribution, this density matrix and its statistical 
sum can be rendered the form of the Euclidean quantum gravity path integral [6] . Its calculation then 
shows if it is dominated either by the contribution of a pure state or a mixed statistical ensemble. 
Thus the dilemma of pure vs mixed state, rather than being postulated, gets solved at the dynamical 
level according to the matter content of the model. 

For models dominated by heavy massive fields this microcanonical ensemble reduces to the pure 
vacuum no-boundary or tunneling state [5], whereas for massless conformally invariant fields the 
situation becomes even more interesting. In this case of the CFT driven cosmology the microcanonical 
ensemble incorporates a possible solution of the cosmological constant problem - the restriction of the 
range of the primordial A by a new quantum gravity scale which is encoded in the conformal anomaly 
of the underlying CFT [5l [6]. Moreover, it contains a mechanism of formation of the red tilted 
CMB spectrum complementary (or maybe even alternative) to the conventional mechanism based on 
primordial vacuum fluctuations in the early inflationary Universe [9] . As it was first observed in [7] this 
follows from a simple fact that thermal corrections to the CMB spectrum enhance its infrared part. 
In connection with this. Professor Dowker might perhaps be interested to know how his pioneering 
results on conformal anomaly [10 and thermal quantum effects |llj in gravity theory underlie the 
foundations of quantum cosmology and, perhaps, have explicit manifestation in the most fundamental 
achievement of contemporary theoretical and observational cosmology - explanation of the large scale 
structure of the Universe. 

In this statistical theory context the operator (jl.ip arises in the one-loop approximation for the 
cosmological statistical sum with t playing the role of the Euclidean time, and the properties of this 
operator essentially differ from those of the Lorentzian dynamics. In the latter case the function g is 
a monotonic function of time because of the monotonically growing cosmological scale factor, whereas 
in the Euclidean case g{T) is periodic just as the scale factor a(r) itself and, moreover, has zeroes at 
turning points of the Euclidean evolution with a = 0, because g{T) oc d(r). This does not lead to a 
singular behavior of F because g also vanishes at the zeroes of g [3] , and the potential term of (II. 1|) 
remains analytic (both g and g simultaneously have a first- order zero). Nevertheless, the calculation 
of various quantities associated with this operator becomes cumbersome due to the roots of gir) - in 
particular, the basis function (|1.3|) becomes singular at each of these roots and cannot be extended 
beyond any of them. Among such quantities is the functional determinant of F which determines the 
one-loop contribution to the statistical sum of the CFT driven cosmology of [3]. Since this operator 
has an obvious zero mode which is the function g{T) itself, the functional determinant of F should, of 
course, be understood as calculated on the subspace of its nonzero modes. The focus of this paper will 
be a collection of issues associated with the calculation of such a restricted functional determinant. It 
will be denoted below by Det*i^. 

Our work largely has a methodological nature and its results in this or that context could be 
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found in numerous papers on mathematical physics. However, their collection below portrays a rather 
illuminating complementarity of various methods which, on the one hand, are focused on calculation of 
this quantity and, on the other hand, embrace rather different fields of mathematical physics, ranging 
from quantization of gauge theories to monodromy method in integrable systems, heat kernel theory, 
spectral geometry, etc. 

The calculation of Det*i^ begins with the remark that there exist several different methods for 
restricted functional determinants. When the whole spectrum of the operator is known this is just 
the product of all non-zero eigenvalues. With the knowledge of only the zero mode, one can use 
the regularization technique [12] or contour integration method [131 Ull [H] to extract the regulated 
zero-mode eigenvalue from the determinant and subsequently take the regularization off. Here we 
use another approach to the definition of Det*F based on the Faddeev-Popov gauge-fixing procedure 
for the path integral in quantization of gauge theories [TB]. In Sect. 2 we interpret the zero mode of 

as a generator of the gauge invariance transformation of the relevant action, so that the reduced 
functional determinant arises as a result gauge-fixed Gaussian path integration. This allows one to 
express it in terms of the Green's function of the operator F on the subspace of its non-degeneracy. 
Remarkably, this Green's function follows from very simple and clear identical transformations under 
the path integration sign, rather than from verbose explanations one usually encounters in numerous 
works on the treatment of soliton or instanton zero modes. 

In Sect. 3 we go over to the calculation of Det*i^ by the monodromy method. First we introduce the 
monodromy of the second basis function ipir) of i^, which is linearly independent of gir). In particular, 
we derive the answer for this monodromy for the so-called multiple nodes case, when the function ^(t) 
within its period range has an arbitrarily high even number 2k of roots, g{Ti) ~ 0, i = 1,2, ...,2k\^ 
The monodromy is presented as an additive sum of contributions of segments of the time variable, 
Ti-i < T < Ti, separating various pairs of neighboring roots. Each contribution is given by a closed 
integral expression in terms of ^(t) on an underlying segment. Then, by the variational method for 
the functional determinant, we express Det*i^ in terms of this monodromy. This actually reproduces 
the known monodromy formula of McKane and Tarlie for the restricted functional determinant with 
periodic boundary conditions [121 1131 HI] , but leaves undefined its overall normalization coefficient 
which is functionally independent of g^r), but can be a function of T - the period of the time range. 

Sect. 4 is devoted to the calculation of this coefficient, and this is a place where various methods 
- WKB approximation, zeta-function technique, heat kernel method, spectral geometry and effect 
of boundaries get efficiently intertwined and complement each other. We begin with the case of the 
nondegenerate operator (|1.1|) subject to Dirichlet boundary conditions on a finite segment of time 
[T_,r_|_], T+ — r_ = T (with a rootless function gir)). It gives the preexponential factor of the time 
evolution operator between r_ and t+, whose dependence on T can be obtained by the combination 
of the C(s)-function method [17] and conformal rescaling incorporating a particular value of the zeta- 
function, C(0), responsible for the conformal anomaly of the theory [10]. The latter is obtained from 
the spectral geometry method and the proper time expansion of the heat kernel [TH [El [20l EU [22] 
in the one-dimensional spacetime with two boundaries at t±. On the other hand, the prefactor of 
this evolution operator with the time-dependent coefficient inclusive is given by the Pauli-van Vleck- 
Morette formula [23^ for the unitary evolution operator (or its Euclidean version [24j). Comparison 
of these two results gives a correct answer for the normalization coefficient. Professor Dowker might 
perhaps be amused to see how this simple effect of boundary terms in the heat kernel expansion 
recovers a correct time dependence of the evolution operator dictated by the Schrodinger equation, 
the situation immeasurably more simple than the one he considered in his pioneering work on boundary 
effects of quantum polarization in curved spacetime [IS] [55]. We accomplish Sect. 4 by applying the 

^ Since a periodic function has only an even number of roots within its period, we will call the case of their lowest 
nonvanishing number, 2k = 2, the single-node one. This is the case of the CFT driven cosmology whose statistical 
sum as a function of the primordial cosmological constant is dominated by the countable set of instantons having k 
oscillations, k = 1,2,..., of the cosmological scale factor a{T) during the Euclidean time period (5] [6] — the so-called 
garlands which carry the multi-node zero mode ^(t) oc a. 
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same zeta-function method to our multi-node case of the operator (jl.ll) subject to periodic boundary 
conditions and prove the T-independent monodromy algorithm for its reduced functional determinant. 
Sect. 5 contains conclusions and the discussion of possible applications of the above formalism. 



2. Zero modes and Faddeev-Popov gauge fixing for reduced 
functional determinants 

When the operator F is nodegenerate its functional determinant can of course be determined by the 
Gaussian functional integral 

(Detx,F)-i/2 ^ poj^g^ y, J Dip cxp^^- S[if]Y (2.1) 
with the quadratic action 



SW]^i^ jjT^{T)F^{T) ^i^jjr l^^' + l'^'j : (2-2) 

where the domain of integration T) and boundary conditions on integration variable ip{t) - one- 
dimensional field - are determined by the class of fields on which the determinant Detpi^ is defined 
(in what follows the label T) will signify both the domain of r and relevant boundary conditions for 
which the operator F is Hermitian). 

The zero mode g{T) of provided it satisfies these boundary conditions, arises as the generator 
of the global gauge invariance of the action (|2.2p under the transformation with a constant e, 



(5^(p = i?(T) e, (2.3) 
i?(r) = ^, \\g\?= I drg^r). (2.4) 

Therefore, for a degenerate operator with the single zero mode the path integral representation of 
its functional determinant ()2.ip can be handled by means of the well-known Faddeev-Popov gauge 
fixing procedure [16]. It consists of imposing the gauge x[v] = and inserting in the path integral 
the relevant Faddeev-Popov factor. This gauge condition xW] and the Faddeev-Popov ghost factor 
Q/llgll can be chosen in the form 

XW]^ I dTk{T)ip{T), <5-x=M^£, (2.5) 
Jv llffll 

Q= [ dTk{T)g{T), (2.6) 
Jv 

where ^(t) is a gauge fixing function and the generator (|2.4[) is normalized to unity with respect to 
inner product on "D. Thus, integration over ip takes the form of the Gaussian functional integral 
with the delta-function type gauge 

(Det;F)-i/2 = const X J Dip s(^x[^]) ^^cxp { - Si^]}, (2.7) 

and serves as the definition of the restricted functional determinant of F. This definition is in fact 
independent of the choice of gauge by the usual gauge independence mechanism for the Faddeev-Popov 
integral. In particular, enforcing the gauge x = means that the field ip is functionally orthogonal to 
the gauge fixing function A:(t) in the metric on V, and the above definition is independent of the 
choice of this gauge fixing function. 

The normalization of the generator (|2.4I) has the following explanation. In local gauge theories the 
Faddeev-Popov path integral is not invariant under arbitrary rescalings of gauge generators. Their 
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normalization is always implicitly fixed by the requirement of locality and the unit coefficient of the 
time-derivative term in the gauge transformation of Lagrange multipliers, Re ^ 1 x e + (which 
follows from the canonical quantization underlying the Hamiltonian version of the Faddeev-Popov 
path integral) . For the global symmetry of (|2.2[) we do not have a counterpart in canonical formalism, 
and such a founding principle as canonical quantization does not seem to be available. Therefore, we 
choose this normalization with respect to unit norm corresponding to the canonical normalization 
of the variable ip in ()2.2p . From the viewpoint of the definition of Det^i^ as the product of operator 
eigenvalues, this corresponds to the omission of a zero eigenvalue of F, 



Det^F 



Ha, 



(2.8) 
(2.9) 



This follows from the orthogonal decomposition of the integration variable ^p{t) in the series of eigen- 
functions >f\{T) satisfying 



dripxir) (Px'{t) = Sxy 



(2.10) 



Representing the delta function of the gauge condition in (j2.7p via the integral over the Lagrangian 
multiplier tt we get the Gaussian path integral over the periodic function 95 (r) and the numerical 
variable tt, 



(Det^F) 



"^/^ = const X Q 

= const X Q \\g\\-'^(Det¥^ 



Dip (It: exp 

-1/2 



S'cff['^(T),7r" 



(2.11) 



Here S'off [(p(r), tt] is the effective action of these variables and F is the matrix valued Hessian of this 
action with respect to = (<p(t), tt). 



Scs['p{t),tt] ~ S[(p] — in / drkip, 

Jv 



F 



6<P S<P' 



F5{t,t') ~ik{T) 

-ik{T') 



(2.12) 
(2.13) 



(note the position of time entries associated with the variables <P = {(p{T),Tr) and <P' = ((/3(r'), tt)). 

The dependence of this determinant on g^r) and ^(t) can be found from its variation with respect 
to these functions. From (|2.1ip we have 



S In (j}et^ F 



-251nQ + 2(51n||.g|| + Tr(jF 



where G is the Green's function of F, F( 
block-structure of (|2.13p is defined as 



(2.14) 



= I and the functional trace of any matrix with the 



Tr 



A{t,t') B[t) 
B{t') a 



dr A{t, t) -f a . 



•D 



The block structure of the matrix Green's function G has the form 



G{t,t') 



Q 




(2.15) 
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where the Green's function G{t,t') in the diagonal block satisfies the system of equations 

FGir,r')=Sir,r')-^^^^, (2.16) 



dTg{T)G{T,T') = 0, (2.17) 

/V 

which uniquely fix it. The second equation imposes the needed gauge, whereas the right hand side of 
the first equation implies that G(r, r') is the inverse of the operator F on the subspace orthogonal to 
its zero mode. 

The trace of the functional block-structure matrix in (|2.14|) corresponding to the variation of g{T) 
reads 

Tr((5gFG) =Tr((5FG(r,r')) = J dT5FG{T,T') _ . (2.18) 

A similar variation of the gauge-fixing function gives a vanishing answer 8k In (Det|, F) = as, of 
course, it should be in view of the gauge independent nature of the Faddeev-Popov path integralU 
This guarantees the uniqueness of the definition of the reduced determinant Det^ F . 

3. Periodic boundary conditions and the monodromy method 

Here we consider periodic boundary conditions for the operator (II. ip which is defined on a circle range 
of the time variable I? = 5^ having the circumferance T . It is parameterized by r 

ro<T<To+r, (3.1) 

with the points tq and tq + T being identified, so that integration over this range will be denoted by 

dT(...) = (fdr {...). (3.2) 



■D 



This range can be infinitely extended to the whole axis — oo < r < oo, multiple covering of 5^, on 
which the function ^(t) and, consequently, the operator F are periodic with the period T, 

g{T + T)=g{r). (3.3) 

The problem of major interest here will be the so-called multi-node case, motivated as it was 
mentioned in Introduction by applications in cosmology, when the periodic function g(r) is oscillating 
and has within its period 2fc simple roots 

To < Tl < r2 < ...T2fe = To + T, (3.4) 

5(n)=0, <?(t,)^0, (3.5) 
gin) = 0. (3.6) 

For simplicity we assume that one of them coincides with the final (or starting) point of this period. 
Another important assumption is that the second order derivative of this function at its roots is 
vanishing, which will be important for analyticity properties of our formalism. 

Another important property of the operator F is its Wronskian relation. For any two functions ipi 
and Lp2 this operator determines their Wronskian ip2] = ^pi^P2 ~ ^\'^2 which enters the relation 

J dT(pi'fif2^ J dTipi^ip2-W[Lpi,ip2] ^. (3.7) 



^In the works involving the treatment of soliton and instanton zero modes it is implicitly assumed that the gauge- 
fixing function coincides with the zero mode itself, A:{t) = (/(t), which considerably simplifies the formalism, but makes 
it less flexible. 



6 



Arrows here denote the direction of action of the operator i. e. (piF — (Fipi), and the Wronskians 
appear as total derivative terms generated by integration by parts of the derivatives in F. When 
both ipi and ip2 satisfy a homogeneous equation with the operator F, their Wronskian turns out to 
be constant. Also the vanishing Wronskian implies linear dependence of these solutions. 



3.1. Monodromy for the multi-node case 

Let us now consider the solution of the homogeneous equation iP{t) normalized by a unit value of its 
Wronskian with g 

F^(t)=0, W[g,i;] = l. (3.8) 

Together with ^(t) this solution forms a set of linearly independent basis functions of F. However, in 
contrast to ^(t) the basis function is not periodic, because we assume that the operator (|l.ip has only 
one periodic zero mode smoothly defined on a circle p. II) . On the other hand, when considered on 
the full axis of r, due to periodicity of g{T) this operator is also periodic F{t + T) ~ F{t). Therefore 
V'(t + T) is also a solution of the equation F{t)^{t + T) = 0, and consequently it can be decomposed 
into a linear combination of the original two basis functions with constant coefficients 

V'(r + T)^V'(T)+A.g(T). (3.9) 

The unit coefficient in the first term follows from the conservation in time of the Wronskian of any 
two solutions of the equation p. 81) . periodicity of (/(r) and an obvious fact that Wlg^g] — and 

1 = W[g{T + T),^{t + T)] = W[g{T), xI:{t) + A5(t)]. (3.10) 

The coefficient A in the second term of p.9|) is nontrivial - this is the monodromy of ?A(t) which will 
play a central role in the construction of the determinant. 

The function iI}{t) can be composed of the set of functions ^iir) defined by ()1.3p on various 
segments of r-range connecting the pairs of neighboring roots of 5(Tjf| 



Mr)=9{r)f T,_i <T,r* <T, , z = l,...,2fc (3.11) 

Jr- 9 iy) 

Here t* are the auxiliary points arbitrarily chosen in the same segments, and all these solutions are 
normalized by the unit Wronskian with gir), W[g,^i\ = 1. The main property of these functions 
S'i(r) is that each of them is defined in the i-th segment of the full period of t where the integral 
(|3.1ip is convergent because the roots of gir) do not occur in the integration range. Its limits are well 
defined also at the boundaries of this segment, 

mn^i) = -^7^, ^'.(T,)--T7^, (3.12) 

because the factor ^(r) tending to zero compensates for the divergence of the integral at r — > — 
and T — >■ Ti-i + 0. 

For an arbitrary choice of auxiliary points r* in (|3.1ip the composite function 

V'(r) =<Z',(r), r,_i<r<r„ (3.13) 

will be continuous in view of (|3.12|) . but the continuity of its derivative will generally be broken, 
because generally the equality 'I'i{Ti — 0) = !Z'i+i(Ti + 0) is not satisfied. However, this equality for 
i = l,2,...,2fc — 1 can be enforced by a special choice of these auxiliary points r*, becoming the 
equation for their determination. The solution for t* is unique, always exists and belongs to the 



^For the extended range 113.191 1 the missing !?'o(t) can be defined by identifying t— i with T2k—i — T and choosing 
some Tq in t_i < Tq < tq. 
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corresponding segment Ti_i < r* < r^lfl On the other hand, the continuity of the derivative of ipiT) 
cannot be attained at ah roots of ^(t), i = 1,2, ...2k, because it would correspond to the existence of 
the second zero mode periodic on the circle, which is ruled out by construction. Therefore, the second 
basis function of F is not periodic on the circle, but in view of periodicity of the operator it satisfies 
the fundamental monodromy property p.9|) . In the next subsection we construct the periodic Green's 
function of F in terms of this monodromy parameter A, whereas here we give in a closed form the 
analytic expression for A as a functional of gir). 

From the definition of the monodromy parameter (j3.9p it follows that in the limit r — > tq , r + T — )■ 

T2fc, 

5(to) \g{T2k) g[T2k-i) J \9[n) g(TQ) J 

where we took into account that g{T2k) = 9{to)- Then the monodromy reads as the additive sum of 
contributions of pairs of neighboring roots of ^(t) [27], 

2k 

A = ^A, (3.15) 
A, = - ^^^^ = ~U,{n)Mn)~^l^,{n^^)^,{n^S]- (3.16) 



9{n) gin-i) 

Because of (13. 6p the functions ^'^(t) are differentiable in these limits, and all the quantities which 
enter the algorithm p.l6p are well defined. In particular, for any such time segment [ri^i, Ti] = [r_, r+] 
the derivatives of ^(t) at its boundaries are given by the convergent integral 

r * 

Note that the integrand here is finite at i/ — r± because of g{T±) — 0. These properties of '^^(t) 
guarantee that the obtained result is independent of the choice of the auxiliary point r* for each A^, 
and the monodromy (jS.lSp is uniquely defined. 

It is important that unlike in the construction of the function ^p{T) which has to be smooth on at 
all roots Ti except tq (the property that was attained above by a special choice of the auxiliary points 
T*), the derivatives of neighboring functions !f'i(T) in p.l5|) - ()3.16|) should not necessarily be matched 
at these junction points. This is because the partial contributions A^ to the overall monodromy A 
are individually independent of r* , 

^^0. (3.8, 

which can be easily verified by using a simple relation (M'i{T)/dT* — —g{T)/g'^{T*). Thus, the mon- 
odromy is uniquely defined and independent of the choice of the auxiliary points r* necessarily entering 
the definition of functions ^'^(t) in Eq. (|3.1ip . 



3.2. Periodic Green's function and the variation of the determinant 

For the calculation of the variation ()2.18p above we need the Green's function of the problem ()2.16|) - 
(|2.17p which should be periodic on the circle p.l|) . To achieve this property we will slightly extend 

^Indeed, the quantity d'j/i(Ti) / dr* = ~ g{Ti) / {t*) is a sign definite function of t* nowhere vanishing on this 
segment, its absolute value quadratically divergent to oo at its boundaries. This in its turn means that '/'^(r,;) is a 
monotonic function of r* which also ranges between —oo and +oo and therefore guarantees the unique solution for t* 
on this segment. 
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the circle domain to the left of the point tq 

To-e <T <Ta+T, e > 0, (3.19) 

with an arbitrarily small positive e and demand that the monodromy of G(t, t') is vanishing for this 
small £-range of r near tq 

G{t + T,t')-G{t,t')=0, Ta-e<T<To. (3.20) 

The ansatz for G{t, t') can be as usual built with the aid of two linearly independent basis functions 
of the operator. One basis function coincides with the periodic zero mode (/(t) and another one is given 
by the function ^{t) built above. Thus it can be represented as a sum of the particular solution of the 
inhomogeneous equation (|2.16|) and the bilinear combination of g{T) and ?A(t) with the coefficients 
providing the periodicity property p.20p . As shown in ^Tj\ it reads 

G(r, r') = Gp(t, t') + i r!(r, r') + a H^^{t, t') + P H^g{T, r') + 7 Hggir, r') , (3.21) 

where 

Gp{r,r')^l{girMT')-^P{r)g{r')) 9{r - r') + ^ (V(r)g(r') - <?(r)^(r')) e{r' - r) , (3.22) 

n{T, t') = w(t) 5(t') + g{T) uj{t') , (3.23) 

H^^(t,t') = ^(r)V(T') , (3.24) 

H^g{T, t') = ^(r) .g(r') + ^(t) V(t') , (3.25) 

Hggiry)^ g{T)g{T'), (3.26) 

and the function uj{t) is defined by 

T T 

lj{t) = ^(r) y dyg{y)kiy) - gir) J dy4'{y)Hy) (3.27) 

with an arbitrary r, . The first term of p.2ip generates the delta- function in the right hand side of 
the equation (|2.16|) . the second term Q,{t,t')/Q gives —k{T) g{T')/Q, while H^^{t,t'), H^g{T^T') and 
Hgg[T,T') represent symmetric solutions of the homogeneous equation with coefficients a, /3 and 7 
which are fixed by the periodicity condition and the gauge condition (|2.17p . 

T,+T 

7-^ J dy^{y)k{y)--^l Jdy^{y)k{y)\ - J dyoj{y)k{y). (3.29) 

The calculation of the variational term p.lSp is based on integration by parts and a systematic 
use of the Wronskian relation p.7p together with equations for g and i/j. The result of this calculation 
is presented in much detail in 27' and reads 



(51n(^Det»Fj = 2 5 In | j^H + (5 In A. (3.30) 
It finally gives the explicit answer for Det, F 

Bet*siF = CiT)xA<l^dTg^{T). (3.31) 
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In fact, this is the McKane-Tarhe formula (Eq.(5.2) of [T^]) obtained by the regularization and 
contour integration methods [121 UM [T5] based on the earher work of Forman [55] for a generic second 
order differential operator. We have reproduced this formula by the variational method for functional 
determinants. Beyond this, the structure of the operator (|l.ip makes the problem exactly solvable 
and gives the monodromy p.9|) in quadratures as an explicit functional of (?(t). 

Final comment concerns the overall normalization in p.3ip . The formula of McKane-Tarlie [T2l 
I14|,ll5j or the variational method which we use below, in principle, give only the ratio of determinants 
for two different operators with different functions g, whereas each determinant contains an infinite 
numerical factor generated by UV divergent product of eigenvalues of F. This factor is independent of 
g{T) but depends on the UV regularization and can be a function of the period T - the only remaining 
free parameter of the problem. In the next section we derive it by the Dowker ^-function method of 



4. (^-function method and zero modes 

4.1. Dirichlet problem and the effect of boundaries 

We begin the consideration of the ^-function method with the case of the functional determinant 
subject to Dirichlet boundary conditions on the initial and final points t± of the time interval of the 
length T 

I?==[T_,r+], T+-T_=r. (4.1) 

Here the focus of our attention will not be zero modes of the operator, but rather - a simple but 
illuminating effect of boundaries, which in the quantum mechanical (that is in one-dimensional) case 
is very universal because it applies to operators of a general form. Thus, we assume that the operator 
(jl.ip is nondegenerate, and its function g{T) nowhere on V equals zero. Therefore, it is not a zero 
mode of the Dirichlet problem because g{T±) = g±^{). 

The Dirichlet problem arises when one considers a semiclassical approximation for the kernel of 
the unitary evolution operator. This kernel, which is given by the path integral, on the one hand 
reduces in the subleading order of this expansion to the one-loop result 

j D.^ exp| - S'[.^] I = const X (^DetuF^ ^^^xpj-S'j, (4.2) 

¥'(t±)=¥'± 

where S is the on-shell action (|2.2p calculated on the solution of classical equations of motion 
(/?(t) = 'p{t\lp±) interpolating between initial and final configurations ip± at t± (the principle Hamilton 
function) 

5^5(r±,^±) = 5[^]|^(^^^^(^l^^) (4.3) 

and Detpi^ is the the functional determinant subject to Dirichlet boundary conditions on [r_,T+]. 

On the other hand, it is given by the Pauli-van Vleck-Morette formula as the semiclassical solution 
of the Schroedinger equation [231 124j^ corresponding to the action (j2.2p 

/ ^^-P{-^[^]}=(-2l^)''''=^H-^}- (4-^) 

¥'(t±)=V± 

This implies the following equality between the functional determinant and the Pauli-van Vleck prcf- 
actor 

DetD-F = const x (4.5) 

\dip+dip-J 



To be more precise, the Euclidean version of this equation corresponding to the imaginary time r. 
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with a constant proportionality coefficient, which is independent of t± because the time evolution 
Schroedinger equation does not admit any flexibility in time-dependent normalization of its solution. 
Below we prove this fact by using the ^-function calculation of the left-hand side of this relation and 
comparing it with the explicit time dependence of the right-hand side - the exercise illuminating the 
role of boundary terms. 

First we find the dependence on time of the expression in the right-hand side of ()4.5|) . For this we 
note that dS / dip+ is the canonical momentum conjugated to the Lagrangian variable ^{t) at t_|_ or 
(/9(r+). This allows one to write a set of relations 



dS , , ^ 



(4.6) 

r— T-f 



u{r+), (4.7) 



(4.) 

where u{t) is a particular basis function of the operator F with the special boundary conditions (which 
follow from obvious boundary conditions for the classical solution with 'f{T^\Lp±) — Lp^^ '^{t-W±) — 

Fu{t) = 0, (4.9) 
m(t+) = 0, m(t_) ^ 1. (4.10) 

This basis function reads 

uir) = 9ir) I 4h f.9(r_) , (4.11) 



so that Eq. (|4.5|) gives 

Detz,[._,.+]F = C(t±) X 5(t+)5(t-) J ^ (4.12) 

where we denoted the normalization coefficient as an unknown function C{t±). In what follows we 
will simplify notations and, without loss of generality, put t_ = 0, = T, C{t±) = C{T) and look 
for the dependence on a single parameter T. 

Now we go over to the zeta-function calculation of the functional determinant. The determinant 
of the nondegenerate operator as a product of its eigenvalues, 

BetvF = YlX, (4.13) 

A 

F^A(r) = A(^a(t), (Pa(t±)=0, reP, (4.14) 

can be expressed in terms of the derivative of the generalized zeta-function which is just the functional 
trace of the inverse s-th power of this operator [T7] 

InDetpF = ^InA = -C'(0), (4.15) 

A 

C(^)=y^i=Tr-L- / , ■ (4-16) 



^ F' J J, F" 



This function is well defined for sufficiently high positive s and exists as analytic continuation in the 
vicinity of s = 0. 



11 



Under the rescaling of the operator or its eigenvalues by a constant coefRcient cr the functional 
determinant transforms as 



Detx, (cr^F) = Yl i'^^^) = ^^'^(o^Detp F, 



(4.17) 



where C(O) plays the role of the number of eigenmodes of F regulated by the zeta-function method 
[To] . Now we use this property to derive the T-dependence of Detp F. 

For this purpose introduce the new operator F„ of the form (jl.ll) with the new function ga{T) 
defined on a rescaled time domain Vc , 



r e [0,r/a] =2?,, g,{T) = g{aT). 



dr^ 9<t{t)' 

Since gaiT) = (""'") — '^^gi'^T') this operator reads 



(j'^dr^ g{<JT) 



dy'^ ' giy) 



y = ar 



e [0,T] 



(4.18) 



(4.19) 



and has as a spectrum the set of eigenfunctions ip^^\T) on P^r, satisfying the same Dirichlet boundary 
conditions and related to the original spectrum of (|4.13p - (|4.14p by 



F,^[1\t) = X,^[l\r), rev,. 
Therefore Detv^F, = Ha equals 

Then, applying (I4.12p to the left hand side of this relation 



T/a 



C{T/a)g,{T/<j)g,{Q) 



dy C{T/a) 



gUy) 



9{T)g{0) 



dy 

g'^iy) 



(4.20) 
(4.21) 

(4.22) 
(4.23) 



and doing the same with the right hand side we obtain the equation for the coefficient function C(T) 

C{T/a) = (7i+2?(0)C(r), (4.24) 

which implies that C(T) = const x r-i-2C(o)^ 

To find ^(0) we can use the heat kernel method in the approximation of the inverse mass (or proper 
time) expansion with the regulator mass parameter added to the operator F. For the differential 
operator acting in a c?-dimensional spacetime with coordinates x it reads 



Crn(s) = Tr 



oo 

jrpj J dtt'-'^Tr exp | - i(F + m2)| 



{4:TT)d/ 



^ oo 
■ \d/2 $Z 



rls + 



n — d\ 



ris) 



(4.25) 



where we have used a well-known Schwinger-DeWitt (or Gilkey-Seely) proper-time expansion for the 
heat kernel trace [H [H [20l [211 [22] 



Tr exp|-i(F + TO2)| 



(4.26) 



n=0 
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Here A„ represent local integrals over spacetime domain B and its boundaries b = 9B, which are 
built of coefficients of the differential operator F. Important property of these coefficients is that for 
odd n they are given exclusively by surface terms and, therefore, identically vanishing for problems 
with a compact spacetime without boundary. For the operator of the form 

F = -a + P, □^.g^'^V^V,, (4.27) 

with the covariant d'Alembertian □ acting in the d-dimensional spacetime with coordinates x = 
and a generic potential term P{x), the first two lowest order coefficients in the Dirichlet case are 

B 

A^^-^ j d'^-^x,/^), (4.28) 

b 

where the Riemann integration measures are built with respect to d-dimensional metric and {d — 1)- 
dimensional induced metric of the boundaryl^ For our case of int erest d = 1 the latter reduces to 
a trivial contribution of zero-dimensional points, /b d'^^^x \/ =2x1, where the coefficient 2 
signifies two end points of B = I? = [0, T]. 

From ()4.25p the value of Cm(0)i which is always given by the first d terms of the local Schwinger- 
DeWitt expansion (analytic in m with < d), in the case of d = 1 reduces to the boundary term A\ 
and reads 

C(0) = U(0) =^ = -i. (4.29) 

m— >0 \/47r ^ 

Therefore 1 + 2^(0) — 0, and C(T) is a constant independent of T, which fully confirms the Pauli-van 
Vleck-Morette formula (14.51). 



4.2. Periodic boundary conditions 

For the periodic boundary conditions with the multi-node zero mode g{T) in the domain 

I? = = To < r < T, (4.30) 

= To < ri < ... < r2fc = T, (4.31) 

the reduced functional determinant is determined by Eqs. (|2.8| )- (|2.9[) which, similarly to ()4.15p bring 
us to the use of the reduced zeta-function 

C(-) = E^- (4-32) 

The scaling behavior is again determined by C*(0) and reads 

Detl, (a^F) = [| (a^A) = a^C-Woet^ F, (4.33) 

which allows us to repeat the steps of the previous section. 

Consider again the operator (j4.18p on a rescaled circle with the rescaled set of zeroes of gaij), 

9a{T-) = 0, 

V, = Sl, = To < T < -, (4.34) 

(7 

= ro-<Tf <... <r2-, = -, = -. (4.35) 

a a 



^In the case of Neumann boundary conditions Ai has an opposite sign on relevant boundaries. 



13 



Its reduced functional determinant is given by Eq. p.3ip 

Det^^ = C (^^^ A'^ £ dTgl{T)^^c{^^ ^jjrg^r). (4.36) 

Here we took into account simple expressions for the monodromy and the zero mode norm of the 
rescaled operator 

= (4.37) 
a 

drglir)^- i drcf{T), (4.38) 

which in their turn follow from the following obvious relations for its partial basis functions ^i{T) 
and its monodromy constituents Af in A"' — Af, 



^:{T)^g,{T) / -^ = lMf,[aT), TU<r<T'^, * = l,...,2fc, (4.39) 

rAv) cr 



Af = -wfm 



A, 



(4.40) 



Applying (|4.36p and p.3ip respectively to the left-hand and right-hand sides of Eq. (|4.33p we get 
the equation for C(T) analogous to (|4.24p . 

C(T/cr) = (72+2C* (0)^(2.)^ (4 41) 

For the calculation of C*(0) again consider the generalized zeta-function regulated by a suffi- 
ciently large mass parameter, so that all eigenvalues A -|- rr? are positive, 

= E / 2V- ^ E TT^^ + 4r ^ + (4-42) 

whence 

C(0) = Cm(0) - 1. (4.43) 

Note that Cm(s) here is not the reduced one, and it contains the contribution of A = shifted by a 
large to the positive range. Therefore, the heat kernel representation applies to it, and the value of 
Cm(0) is again given by Eq. (|4.25P . but this time for a compact domain V = S"" without boundaries. In 
this case all odd number coefficients are vanishing A2k+i = 0, because they are exclusively contributed 
by boundary terms, and Cm(0) = 0. Therefore, C*(0) = C„(0) = -1 and the coefficient C'{T) in (|OT|) 
and p.3ip is again a constant independent of T. 

Its actual value within zeta-function regularization can be determined for a particular case of the 
constant function g{T) = c corresponding to the operator F = —d^/dr^ with the explicit spectrum of 
eigenfunctions and respective eigenvalues 

(^0 = 1, Ao = 0, (4.44) 



(pinir) ^sini^—Tj , ip2n{T) ^ cosI^-j^tJ , ^ j , n=l,2,... . (4.45) 

The logarithm of the corresponding restricted determinant - the product of all nonvanishing eigen- 
values regularized by zeta-function method - equals 

lnDet,(^-|l) = 2£ln (^^)' = 41n (^^) UO) ~ 4Cj,,(0) = 21nT. (4.46) 
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Here C-r(s) = X^J^i Riemann zeta function having the following particular value C-r(O) = — 5 

and the value of its derivative Cfl(O) — —^\n2Tr. On the other hand, the basis functions and the 
monodromy A for this operator read 

= T, 

A = ^. (4.47) 

Therefore according to p.3ip Det, (-d^/dr^) = C{T)T'^, and the comparison with (|4.46p gives the 
T- independent result C{T) = 1 for the normalization coefficient in (|3.3ip . 

5. Conclusions 

Thus we see that the above combination of methods gives exhaustive answer for the reduced functional 
determinant of the operator having a multi-node zero mode in the periodic boundary value prob- 
lem. This determinant expresses in terms of the monodromy of its basis function, which is obtained in 
quadratures as a sum of contributions (I3.15l) - (|3.16l) of time segments connecting neighboring pairs of 
the zero mode roots within the period range. Few words are in order here just to reiterate our special 
interest in this particular problem, briefly mentioned in Introduction. 

The operator F determines the one-loop statistical sum for the microcanonical ensemble in cosmol- 
ogy generated by a conformal field theory [SI [SJ [3] . This ensemble realizes the concept of cosmological 
initial conditions by generalizing the notion of the no-boundary wavefunction of the Universe to the 
level of a special quasi-thermal state which is dominated by instantons with an oscillating cosmological 
scale factor a{T) of their Euclidean FRW metric. These oscillations result in the multi-node nature of 
the zero mode g(r) ^ a(r) of i^, which itself arises as the residual conformal Killing symmetry of the 
FRW background. This, in particular, explains the motivation for the gauge-fixing treatment of the 
zero mode considered above. 

As was mentioned above, a very attractive feature of the cosmological microcanonical ensemble 
is that in the case of the CFT driven cosmology it suggests a possible solution of the cosmological 
constant problem - the restriction of the range of the primordial A by a new quantum gravity scale, its 
value being encoded in the conformal anomaly of the underlying CFT [316]. Moreover, as suggested 
in [7], these microcanonical initial conditions admit inflationary scenario in the early Universe and 
can provide a thermal input in the red tilt of the COBE part of the CMB spectrum. This tilt can 
be additional or, perhaps, even alternative to the conventional red tilt generated from primordial 
vacuum fluctuations of 9J . This makes the hypothesis of microcanonical initial conditions in quantum 
cosmology not only feasible, but also obscrvationally verifiable, perhaps, in a foreseeable future. Also, 
the statistical sum of this ensemble is likely to predict interesting phase transitions for multi-node 
cosmological instantons [29] which makes physics of this model very rich and interesting. The results 
and methods presented above seem indispensable for a further progress in these intriguing issues. 
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